A new theoretical approach to 1:1 electrolytes at low temperature 
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A new theoretical approach to 1:1 electrolytes at low temperature is developed, RPM and SAPM 
are studied with this approach, and their critical points of first order phase transition are calculated. 
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This report presents a new theoretical approach to 
ionic systems at low temperature. The systems un- 
der study are primitive models of electrolytes: Equal 
numbers, N, of positively and negatively charged hard 
spheres of diameter a+ and a_ in a volume V, carry- 
ing charges +qo and —qo respectively, interact via the 
Coulomb potential ±g§ /Dr in a medium of dielectric con- 
stant D. The size asymmetry of + and — ions is mea- 
sured by A — a^/a-\-, with A = 1 being the special case of 
the Restricted Primitive model (RPM) and otherwise the 
Size Asymmetric Primitive Model (SAPM). Such systems 
undergo a first order phase transition at low temperature 
and moderately low density, they have received increas- 
ing attention in recent years (see, e.g. In 
the following, we define a = \ {a+ + a_ ) and the normal- 
ized reciprocal temperature (3 = qQ/{kBTD), which has 
a length dimension, we denote the total particle number 
density by p, and because of the symmetry with respect 
to the exchange of -I- and — ions, only A < 1 is considered. 

There are 2iV^ — N pairwise interactions in the sys- 
tem, among which are negative (between unlike ions), 
iV^ — N are positive (between like ions). Intuitively, if 
the N'^ — N positive interactions could be used to can- 
cel N'^ — N negative interactions, we would only have 
N negative pairwise interactions left in the system. And 
this is desirable because it's much more manageable to 
deal with N than 27V^ — N interactions, analytically or 
numerically. In this context, let's review a pairing proce- 
dure first introduced by Stillinger and Lovett 6] . In their 
procedure, all the distances between two unlike ions are 
computed, the first pair is defined as the two unlike ions 
that have the closest distance, and this step is repeated, 
taking into account only ions that remain unpaired, un- 
til all the ions in the system are exhausted. If a certain 
prescription Q is used on the situations (of zero weight) 
when one ion has more than one unlike ion at the same 
distance, the result of such procedure is that each posi- 
tive ion has one and only one negative ion as its partner, 
and vice versa. We call this procedure "Closest Pair- 
ing" (CP). The following numerical study was carried 
out to exam the energy profile of CP pairs: N — 500, 
a+ = a_ = 1, ions are randomly |3| placed in a cubic 
box of volume V, with hard wall boundary condition; 
pairs are formed through CP procedure, and the N at- 
tractive inter-ion interactions of CP pairs are summed 



at p = 0.05. cross: with no restrictions on the selection of 
configurations; circle: with DNR imposed on the selection of 
configurations; dashed line: /, as reference. 



up, divided by N, denoted by Icp, and then compared 
with /: the sum of all the 2N'^ — N pairwise interactions 
divided by N. We found that for a wide range of densi- 
ties, at least for p e [0.01,0.2], Icp is within close range 
of / for configurations with / e [— 1,— p^/'^]; while for 



Icp has a value less but around —p" 



1/3 and 

shows no significant change for different configurations; 
and of course, Icp cannot have values below (— 1/a), here, 
-1. A typical plot of Icp vs. I is shown in Figure ^ 

It is then adequate to approximate / by Icp in the low 
temperature regime: (3 ^ a, when, expectedly, most un- 
like ions are closely paired up and the value of / is in 
the close neighborhood of (—1/a). For clarity of pre- 
sentation, we call such approximation "CP Approxima- 
tion (CPA)". But configurations with / > -p^/^ 
pose a potential problem for us. Although in the ther- 
modynamic limit when N,V — > oo, such configurations 
would have virtually no contribution to any physical phe- 
nomenon we are interested here, they do occupy a rela- 
tively very large volume in phase space because they al- 
low particles to move much more freely. Thus if no pre- 
cautions are taken in the calculation of the partition func- 
tion with / approximated by Icp, this large but ought to 
be insignificant phase space volume would be associated 
with a significantly magnified Boltzmann factor, and this 
could significantly affect the result quantitatively, if not 
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qualitatively. In order to prevent this from happening, 
we put a restriction on the selection of the allowed con- 
figurations, called "Detailed Negativeness Requirement" 
(DNR): two like ions can not get closer than a distance 
S) unless at least one of them has at least one unlike ion 
at a position that is closer than S). A similar numerical 
study was carried out in which only the configurations 
that comply with DNR are sampled, and we found that 
DNR does eliminate unwanted configurations 0; a cor- 
responding plot of Icp vs. I is also presented in Figure 
(we see that DNR also limits local energy fluctuations 
and admittedly deserves further study). 

In the following, we present an analytical scheme that 
employs CPA with DNR. 

In the Grand Canonical Ensemble, the system is de- 
scribed by five parameters (a. A, z+, z_, /3). The partition 
function is: 



JV+=0 Af_=0 



2 — r. 



^ + V';ic(r-ij,(Tj,(Tj|a, A)]} 



(1) 



where Oi is the sign of ion i and j ~ jr^ — x.j \ is the dis- 
tance between ion i and ion j. The hard core interaction 
V/ic has the usual meaning, here we have to specify the 
signs of the two ions involved because of the size asym- 
metry. 

We now carry out the CP procedure on any configura- 
tion of the system, and group pairs that have the inter- 
ion distance in [a, a-|-^], ^ being infinitesimal; denote the 
number of such pairs by iVg. For the rest of the ions in the 
system, this operation and DNR would change the first 
two system parameters from (a. A) to (a + ^, A), where 
A = min(l, A -I- ^(I -I- A)/a). Now we have: 
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exp{-/3E[ E v,^™ + E^]} 

m— 1 n—1 



(2) 



where X = R © s is the six-component vector comprising 
the position of the center of a pair, R, and the relative 



displacement vector, s, which points from — ion to + ion 
of this pair. In the integration range of X, v is the range 
of s delimited by the requirement that s G [a, a + 
ypsj.m denotes the interaction between pair I and single 
ion TO, written explicitly: 

^ps, 1,771 — ^^s(X;,r^, 

am\a, A) = Vps(R/, s/, r,„, (Tm|a, A) 

= Vhci\i'Rl + |) - rml, +1, A)) 

+ 14e(|(R; - |) -r„|,-l,a„>,A)) 

where Vd is the extra spacial exclusion imposed by CP 
and DNR: 

Krf(R/,s,,r„,a™)-|Q otherwise. 

Notice that in the expression of Vps there are no electric 
interaction terms because of CPA, and this makes Vps a 
short-ranged interaction. Vpp^ „ is the interactions be- 
tween two CP pairs, I and n, which is also short-ranged 
because of CPA. Since ^ is infinitesimal, configurations 
with more than one chosen pair in any finite region can 
be neglected, so can Vpp, and the summation over Nq in 
Equation (0) can be replaced by an exponential form: 



exp{^z_l_z_e" x M} 



(3) 



where 



N++N. 



M = j dH Jdujexp{-I3 ^ yps(R, a(tj), r„, cr„|a. A)} 

V 47r ™=1 

(4) 

Lu is the orientation of the pair supplying the direction of 
a.{uj) whose magnitude is a. 

To evaluate Equation I^J, we decompose the exponen- 
tial form in the usual way: 

N++N- 



M^JdnJdu Y[ (/„, + ! 
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N++N- 



(5) 



m— 1 



n,7i —1 



where: 



/,„ = cxp{-/3yps(R,a(u;),r„,cr„|a, A)} - 1; 
M = JdR JdLul = inV; 

V 4tt 

S,n = JdR JdLufjn; 
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Higher order terms which involve more than two f's can 
be defined similarly. Notice that S is the same for all ions 
of the same sign, in the following we'll use 5+ (a, A) and 
<S_(a, A) to denote them respectively. I?„^„' depends on 
the value of (a, A) and the relative position and signs of 
the ion n and n'; it acts like an extra 2-body interaction 
for the rest of the ions in the system. Similarly higher 
order terms act like extra multi-body interactions. Now 
the partition function becomes: 



S = exp{^z+z_,e» x A/"} 



EE? 



z.z 



N+=0 JV_=0 



xexp{-- + V^.c.(CTi,crj,rij|a + ^,A)]} 

1,3 = 1 

X exp{^2;+z_e" x ^ Vn^' + . . •} 



(6) 



where z± = z± x exp{^z^z_e^ a'^ x iS-|-(a, A)}. 

We see an interesting pattern in Equation the 
LHS denotes the partition function with system param- 
eter (a, A, z+, z_, while the RHS has an exponential 
term followed by the expression of a partition function 
with system parameter (a -I- ^, A, z+, z_, /3) and the ex- 
tra interactions introduced by V and high order terms of 
Ai. We can of course assume that there are also these 
extra interactions on the LHS, but they have the spe- 
cial initial value of 0. Because of this similarity, the 
same procedure we performed on the original partition 
function can be carried out on the new partition func- 
tion on the RHS and then repeated: a differential proce- 
dure that resembles the Renormalization Group Theory, 
except that no fixed point of the system parameters is 
expected. Since Vps is short-ranged, we will neglect V 
and all the higher order terms of A4 in the following 
for simplicity; their effects will be discussed later in this 
report. The system parameters are now interpreted as 
variables to be "scaled up" in this procedure, with no 
confusion, denoted as (a. A, 5+, z_, /?) that take the ini- 
tial values of (a. A, z+,z-,/3), and a serves as the scaling 
factor. From the above derivation, the differential equa- 
tions of the other system parameters with respect to d 
are: 

3^,.^, r (l + A)/a ifA<l; 

[ otherwise; 

|3 (7) 
z'±{d) = z± X z+z^e'Sa? x S±(a, A); 

/3'(a) = 0. 

which are easily solved. This procedure can be carried 
out until the temperature is not considered low anymore 



when a = K (3, and we have: 



S = exp{y / z_|_(a)z_ (a)e ° 47ra da}x 

a 

S(A,A(A),f+(A),l_(A),/5). 



(8) 



where S(A, A(A), z+(A), z_(A), /?) is the contribution 
from "free moving" ions, which can be evaluated by lin- 
earized theories such as the Mean Spherical Approxima- 
tion (MSA) or the Debye-Huckle (DH) Theory. The ex- 
act value of A is, of course, not important; in this report, 
we set A = /3/2. 

Equation (|HJ| concludes the above derivation, with sim- 
plifications of course: recall that 2? and higher order 
terms of have been neglected. If these terms are in- 
cluded in the derivation, theoretically speaking, a simi- 
lar differential procedure can still be constructed, but it 
will have greater complication. Such simplification over- 
estimates the repulsion between the CP pair and sin- 
gle ions at each step of the derivation, consequently, the 
value of pressure is raised; and the contribution from 
S(A, A(A), z+(A), z_ (A), /?) is lower estimated because 
z±(a) is decreasing faster than it ought to be. But if 
the density of the system is not too high, it should not 
significantly affect the result. And if we are interested 
in situations with very low temperature: (3 ^ P~^^^^ 
how it is when the system is close to the critical point of 
the phase transition, the contribution from " free moving" 
ions can be neglected, because, not surprisingly, almost 
all ions would be closely paired up at such a low tem- 
perature; this can also be seen from the fact that, when 
P :§> a and z± is sufficiently big (so that p~^^^ ^ if 
S(A, A(A), z+(A), z_ (A), /3) is evaluated in the same fash- 
ion as above (the integration with respect to d now goes 
from A to cxd), with CPA but without DNR: a clear over- 
estimation, its contribution is still negligible 0- The 
partition function then takes a very simple form, written 
explicitly: 



where z 



S = exp{l// 



and: 



ze s Aira 



1 - z!ag{x)dx 



da} 



gix) = a;'^e^/^[5+(2;,A(2:)) +5_(x,A(a;))] 



(9) 



(10) 



The description of the system by Equation lO, how- 
ever, does not exhibit a phase transition. Physically, we 
have seen that when CP pairs are evenly (randomly) dis- 
tributed, I Pi Icp and the electric interaction between CP 
pairs, denoted by Vppe in the following, overall doesn't 
have significant effects due to cancellation. But this 
ceases to be true for CP pairs in clusters, which are shown 
by MC simulations (e.g. [lli]) to form abundantly when 
the system is close to phase transition. Recall the well- 
known relationship: U = -g^ J E^(r)dr, for two CP pairs, 
lowering Vppe would lower U, which in turn lowers their 
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electric field as a whole and also makes the "shape" of 
their E(r) more concentrated. Because their electric field 
serves as the agent of the interaction between the two 
CP pairs and the rest of the system, lowering Vppe would 
then make them more isolated from the rest of the sys- 
tem and their negative interaction less likely to be can- 
celled by positive interactions. For larger clusters, this 
effect becomes even more significant. In the following, 
we construct a VdW-like term to estimate the effect of 
this " extra" attraction between CP pairs and extend the 
scheme to exhibit phase transition. 

The previous derivation is treated as the leading order 
calculation, from which the mean inter-ion separation of 
CP pairs is found to be (2|: 

= ^ r -"[^-^g(^)(^-i)-^^]^-^ a (11) 
pJa [1 - z J: g{x)dx]^ 

The change of pressure as the result of including this 
"extra" pair attraction is then estimated to be: 

AP2 = ii?2,e((.),A,/?)(^f (12) 

where 



i?2,e((s),A,/3) = 




where Vppe = Vppe{'Ri,uj,'R2,u!2\{s), \) gives the elec- 
tric interaction between two identical CP pairs of inter- 
ion separation (s) for configurations allowed by CP and 
DNR, and otherwise 0. Here, i?2,e accounts for the "ex- 
tra" attraction introduced by 4-ion (two CP pairs) clus- 
ters, and, following the spirit of the above discussion, 
we propose the use of energy standard in defining clus- 
ters: £2 is the maximum electric interaction between two 
pairs for them to be considered as a 4-ion cluster, big- 
ger clusters can be defined in a similar fashion. The val- 
ues of £2 and £s for clusters involving more than two 
CP pairs are chosen so that, in the sense of statisti- 
cal average, when I < (— 1/a) (recall Icp > (— 1/a)), 
[I cp + 1 cluster) " closcly" approximate /, where Iduster de- 
notes the "extra" attraction introduced by clusters. £s 
would of course depend on the temperature and the den- 
sity of the system. Here, we take a short cut and choose 
£2 = — 0.48/(s) so that the critical temperature of RPM 
matches the result of MC simulations, which has con- 
verged to ~ 0.050{qQ /kBDa) in recent years. 

Since the big asymmetry cases raise new problems with 
the formation of special micro-structures |j] and is still 



under investigation, here we focus on situations with only 
moderate asymmetry. Some of the critical points ob- 
tained from this work are reported in Table |ll in which 
T* — ksTcDa/qQ and p* = PcO? are the reduced criti- 
cal temperature and reduced critical density respectively, 
values from the MC simulation of , if available, are 
also listed in the table as comparison, they are denoted 
by T* Y and p* y. 



TABLE I: locations of critical points 



A 


T* 
^ c 


P*c 






1.00 


0.049(6) 


0.070(3) 


0.0492(2) 


0.073(2) 


0.75 


0.049(2) 


0.063(4) 


0.0488(2) 


0.072(2) 


0.50 


0.048(2) 


0.053(5) 


0.0475(3) 


0.070(2) 


0.45 


0.047(7) 


0.051(2) 






0.40 


0.047(0) 


0.048(6) 






0.35 


0.045(9) 


0.046(1) 






0.30 


0.044(3) 


0.043(5) 






0.25 


0.041(6) 


0.040(6) 


0.0422(3) 


0.059(3) 


0.20 


0.036(7) 


0.038(0) 


0.0386(4) 


0.051(3) 



For the critical density of RPM, recent MC simulations 
have converged to p* ~ 0.075, the result of this work (the 
case with A = 1.00) has a discrepancy of less than 10%. 
For the critical temperatures of all the values of A, we 
see a very good agreement (less than 10% of discrepancy) 
with the simulation results of |J| . More importantly, the 
result shows that both critical temperature and critical 
density decrease as the size asymmetry increases, and the 
decrease is small until A ^ 0.4, after which they decrease 
more dramatically, these features are consistent with the 
results of MC simulations 0, 0|; to our knowledge, it has 
not been shown before theoretically. 
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